arXiv:1508.02821v3 [math.DG] 6 Apr 2016 


HARNACK ESTIMATE FOR MEAN CURVATURE FLOW ON 

THE SPHERE 

PAUL BRYAN AND MOHAMMAD N. IVAKI 


Abstract. We consider the evolution of hypersurfaces on the unit sphere 
gn+i their mean curvature. We prove a differential Harnack inequality for 
any weakly convex solution to the mean curvature flow. As an application, 
by applying an Aleksandrov reflection argument, we classify convex, ancient 
solutions of the mean curvature flow on the sphere. 


We study convex liypersurfaces Mq, n > 2, without boundary, which are smoothly 
embedded in Let Fq : M" —>• be a smooth embedding of Mq. The mean 

curvature flow with initial data Fq is a family of hypersurfaces given by embeddings 
F : M” X [0,T) —^ F{-,0) = Fq{-) that moves in the direction of the unit 

normal vector u with speed equal to the mean curvature H, the trace of the second 
fundamental form A(V, V) over the tangent vectors V, e.q., 

dtF{-,t) = 

Here, H{p, t) is the mean curvature of the hypersurface Mt := F{M'^,t) at the point 
F{p,t) where the unit outer normal vector is Let D denote the Levi-Civita 

connection of 

Theorem (Harnack estimate). Suppose Mt is a smooth, weakly convex solution of 
the mean curvature flow on the time interval [0,T). For t > 0 we have 

TT 

9tH+ — + 2DvH + A{V, V)-nH>0 
for all tangent vectors V. 

Comparing this with Hamilton’s Harnack inequality [5] for the mean curvature 
flow in here we gained a bonus term —nH through the ambient space 

This allows us, by employing a parabolic Aleksandrov reflection argument first 
developed in [6, 7, 8], a rather convenient classification of ancient, smooth, weakly 
convex solutions of the mean curvature flow. 

Theorem (Classification of ancient solutions). The only ancient, weakly convex, 
embedded solutions of the mean curvature flow on the sphere are equators and 
shrinking geodesic spheres. 

Let us note that such a classification has been obtained already in [13, The¬ 
orem 6.1], with a relatively short proof applying the maximum principle to the 
quantity (||A|p — As usual for such arguments, the Codazzi equation 

is employed; therefore, the result pertains to n > 2. The first author and Louie 
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[4] used a completely different argument to obtain the classification for n = 1. 
First, a Harnack inequality is obtained, which ensures that the backwards limit as 
t ^ — oo is an equator. Then the Aleksandrov reflection argument is employed to 
show maximal reflection symmetry, immediately resulting in the classihcation. Our 
work here extends the techniques of [4] from n = 1 to n > 2 formally following the 
same procedure of obtaining the Harnack inequality and then using the Aleksan¬ 
drov reflection to show the maximal reflection symmetry. We believe, that unlike 
the methods in [13], our argument is more broadly applicable: once the Harnack 
inequality is obtained for a curvature flow, the remainder of the argument holds 
with very little modification. This theme will be explored for a range of curvature 
flows similar to those studied in [2] in a forthcoming paper. In [2], the Gauss map 
parametrization and support function are used to great effect, considerably simpli¬ 
fying the computations. A genuine difficulty on the sphere is obtaining a similar 
useful parametrization under which we can perform the calculations. 

Before moving on to the details, a few words about the broader context are 
in order. It is well known that ancient solutions arise as singularity models for 
curvature flows (e.g. [12]); therefore, their study is of great importance. On the 
unit sphere, as opposed to the Euclidean case studied in [13], where classification can 
be obtained only with additional assumptions such as pinching, we find a complete 
classification and strong rigidity. The difference is the compactness and positive 
curvature of the sphere. We might conjecture that for the mean curvature flow in a 
positively curved, compact background, there exists at most one non-trivial ancient 
solution emanating from each closed minimal hypersurface. However, in general, it 
may be too much to hope that such ancient solutions exist. 

Classifications have been obtained for the curve shortening flow in the plane [9] 
and the Ricci flow of surfaces [10]. Both examples include self-similar shrinkers 
(round circles for the former and constant curvature spheres for the latter) and 
a unique Type H ancient solution with curvature blowing up faster than jfj as 
t —>■ — cx) (the Angenent oval [paper clip] in the former case and the Rosenau in 
the latter). Our results here rule out the possibility of such Type H singularities 
and show that in any dimension, the only weakly convex, ancient solutions are 
the “self-similar” shrinkers: here, self-similar means conformally equivalent to the 
self-similar family of shrinking spheres in Euclidean space after realizing the sphere 
as (locally) conformally equivalent to Euclidean space. 

In higher dimensions, classification results for ancient solutions in Euclidean 
space are less complete and a wider range of possibilities may occur. As mentioned 
above, the results of [13] require additional pinching assumptions to deduce that 
such ancient solutions are shrinking spheres. Wang [18] studied translating ancient 
solutions and discovered the existence of non-rotationally symmetric translators 
(though they blow down to spheres or cylinders). Furthermore, there are solutions 
asymptotic to “ovals” in that the center looks like a cylinder, but the ends look like 
“bowls” [1, 14, 19]. These solutions have similarities to the Angenent oval in the 
plane. 

In the sphere, we may consider an ancient solution moving by a rotation to be a 
translating solution, but we quickly realize that no such convex solutions can exist 
because convexity and the containment principle force the flow to lie in a fixed 
hemisphere at all times. Our results also show that no non-self-similar ancient 
solutions exist. The phenomena of becoming more and more oval as f —>■ — oo for 
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convex solutions appears to be ruled out by the compactness of the sphere; such 
a solution would have to touch an equator at the “ends” and “fatten” out in the 
middle, which forces it to touch the equator in the middle and thus everywhere (by 
convexity). 

Lastly, the connection between Harnack inequalities, ancient solutions and soli- 
tons deserves mention. In Euclidean space, the classic work [5] espouses the philos¬ 
ophy that Harnack inequality should be an equality on solitons. We have already 
observed that unlike Euclidean space, no translating convex solitons exist in the 
sphere. The equivalent of homothetic solitons may be taken to be those solutions 
of the mean curvature flow following the integral curves of the conformal ’’posi¬ 
tion” vector field sin(d)Vd, where d denotes the distance to a fixed point (the point 
of collapse). This vector field generates a one-parameter subgroup of conformal 
transformations, under which the mean curvature flow is not invariant, though the 
family of shrinking geodesic spheres is such a soliton and hence these conformal 
solitons are of great interest. The lack of conformal invariance suggests that con¬ 
formal solitons do not satisfy equality in the Harnack inequality, and this is indeed 
the case for the shrinking spheres. Interested readers should consult [15, 16, 17] 
for more details on solitons interpreted as the flow along integral curves of a vector 
field, particularly [3] for details on conformal solitons. 
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1. Preliminaries 

Let g = {gij}, A = {h^}, and Rrriijki denote, in order, the induced metric, the 
second fundamental form, and the curvature tensor of M". The mean curvature of 
M” is the trace of the second fundamental form with respect to g, H = g^^hij. Let 
'g and Rnia^-ye denote, respectively, the metric and the curvature tensor of 
Greek indices run through {0, • • • , n} and Latin indices belong to the set {1, • • • , n}. 

Write u for the outer unit normal to M*. For a fixed time, we choose a local 
orthonormal frame {9o = u, • • • , 9^ = ^^0^, • • • , in We use the following 

standard notation 

Iri = g^^h^m 

|A|2 = g^^g^^h^uhij = 

C = g^^9’''9^"h,khimKj = h\hih\. 

Here, is the inverse matrix of {9ij}- 

The relations between A, i?TO, and Rm are given by the Gaufi and Godazzi 
equations: 

Rfn^jki — RTTtijj^i hikhji huhj]^ 

V ihjk — V kdij . 
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Moreover, Vi and A commute as follows 

(V,A-AV,)/ = 

for all smooth functions on M”. Therefore, in view of Rap-ye = a-yVpe ~9aeVp-y)^ 

A = 1, and the Gaufi equation we have 

(1.1) (V,A - AV,)i7 = {{h^)T - HhT -{n- l)6rX)VmH. 

2. Evolution equations 

In this section, we assume that Mt is a strictly convex solution of the mean 
curvature flow. 

Lemma 1. The following evolution equations hold under the mean curvature flow: 

( 1 ) dtgtj = -2Hhij 

(2) dtg^^ = 2Hg^'^g^'^hij = 2Hh^i 

(3) dthl = ViV,H + + XHSf 

(4) dthi = Ahj + \A\'^hl + X{2H6f - nhfl} 

(5) dthij = Ahij + \A\‘^hij - 2H{h^)ij + X{2Hgij - nhij} 

(6) dtH = AH + H\A\^ + nXH 

(7) (dtA - Adt)H = 2Hh^^y,VjH + 2h^^y,HVjH. 

Proof. The computations for (1 — 6) are straightforward (see for example [11]). To 
obtain (7), we note that 

{dtA - Adt)H = {d^djH - rf^-5fci7) dtg^^ - g^^dkHdtT^^, 

where Tfj = g^^digji + djgu — digij). Since dtP^j is tensorial, using (1) and (2) we 
may carry out our calculations in a normal frame to prove the claim. □ 

Lemma 2. Under the mean curvature flow we have 

dfldtH) =AdtH + AHh^^ViVjH + 2h^^\JiHVjH 
+ (|A|2 + nXfldtH) + 2H^C + 2Ai^^ 

and 

VAH = AV,H + Vfl\A\^H) + ((h^)- - HhT)X7mH + XV,H. 

Proof. Using (3), (6) and (7) in Lemma 2, we calculate 

dfldtH) =dt{AH + H\A\^ + nXH) 

=AdtH + 2Hh^^V,VjH + 2hf^V,HVjH 
+ {\A\^+nX){dtH)pH{dt\A\^) 

=AdtH + 2Hh^^V,VjH + 2h^^V,HVjH 
+ {\A\^ + nX){dtH) + H{2h^^ViVjH + 2HC + 2XH^) 

=AdtH + 4:Hh^^V,VjH + 2h^^V,HVjH 
+ {\Af + nX){dtH) + 2H^C + 2XH^. 


HARNACK ESTIMATE FOR MEAN CURVATURE FLOW ON THE SPHERE 


5 


To obtain the second evolution equation, we use identity (1-1) and part (6) of 
Lemma 2: 

y,dtH =y^{AH + \A\^H + nXH) 

+ nXSJ,H 

+ {{h'^)^-HhT-{n-l)5TX)VmH 

=AV,iL + V,{\A\^H) + ((h2)7* - HhT)VraH + XV,H. 

□ 

In the sequel, {6*-^ } denotes the inverse of the second fundamental form. 

Lemma 3. 

dtib^^ViHVjH) < - b^'^b>^AhmnV,HVjH + 2V^VjHAViH 

+ 2K-^ViHVjH + \ A\^b^^ViHVjH + 2Hb^^V ,\A\^V jH 
+ nXb^^V,HVjH. 

Proof. Observe 

-XP^V'^{2Hgmn - nhmn}ViHVjH + 2XP^ViHVjH < nXP^ViHVjH 

and 

dtP^ = - 

= - P’^V^{Ah,nn + \A\^h,nn - 2H{h^)mn) 

- XP^V^{2Hgmn - nhmn] 

= - P'^V^Ahmu -\A\'^b^^ + 2Hg^^ 

- XP"^y>^{2Hgmn - nhmn}. 

Therefore, using the second part of Lemma 2, we obtain 
dt{b^^V,HV,H) = - P^V^{Ahmn + \A\^hmn " 2H Q?)mn)V ,HV, H 

- XP"^V^{2Hgmn - nhmn}V,HV,H 

+ 2P^ {AV,H + V,{\A\^H) + ((h^)™ - HhT)VmH + XV,H)VjH 
< - P^V^Ahmn^iHVjH + 2P^V,HAV,H 

+ \A\^P^V,HVjH + 2HPiVi\A\^VjH 
+ 26*^ ((h^)™ - HhT)VmHVjH 
+ nXP^V ,HV jH 

= - V'^V'^AhmnViHVjH + 2P^VjHAV,H 

+ 2h^^V,HVjH + \A\'^P^V,HVjH + 2Hb^^Vi\A\'^VjH 
+ nXP^ViHVjH. 

□ 

Using the identities 

Vmb^^ = -P’^V^Vmhki 
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and 

= - h^'^V'^Ah^n 

+ + b^>^V^b‘^}gP‘>VpK,V,hM, 

we compute 

A{b^^V,HVjH) = - b^"^V^V^HVjHAhmn + 2b^^VjHAV,H 
+ 2{b^^b'^^b^^}V^HVjHgP<^VphrsVghki 

- 4g^^bH^^S/^hki\/nV,HVjH + 26*^g™"(V^V,i7V„V,if). 
Thus we have proved: 

Lemma 4. 

dt{b^^V,HVjH) <A{b^^V,HVjH) 

-2{b^^b'^^b>^}V^HVjHgPiVphrsVghki 
+ 4.g^^b^'^V‘ymhklyny^HV,H - 25*^g™"(V™V,7?V„V,i7) 
+ 2h^^\JiHVjH + \ A\%^^\JiHVjH + 2Hb^^V jH 

+ nW^y,HVjH. 

Lemma 5. Define 0 := dfH — D^ViHVjH. Then 

dtO > A0 + + (|^|2 ^ 

H 

mn pq 

+ 2{g^%^P - + 2XHf 

where 

Vmn — ^nH A Hifl ')mn b ^r-hb^si^mn- 


Proof. Note that 

, (0-nAH)^ 

H 


=2 


+ 4\A\'^AH - 
H H 


+ 2\A\*H - ViTfVji? 

+ I 


mn pq 

2{g”^%^P - ^-^}Vmnripg 

=2g^%^PX/m^nHypVgH - 2^^^ 

12 

+ V^VjiJ - A\A\^AH 


- 4g'^%’^PD^Vm^nHVrHVshpg + A—P^V,HV,H 

H 

+ 2H'^C -2\A\^H 

- 2HP^ViHVj\A\‘^ + 4|Ap6*^V,iJVjiJ 

+ 2{b"b'^^V^}V,HVjHgP‘iyphr.syghki - ^ {D^V.HVjHf . 
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Now the claim follows from adding up these last two identities and considering 
Lemmas 2 and 4. □ 


3. HARNACK ESTIMATE AND BACKWARDS CONVERGENCE 

3.1. Harnack estimate. If Mq is not an equator, the strong parabolic maximum 
principle and the evolution equation of in Lemma 2 imply that for any t > 0, 
Mt is strictly convex. For strictly convex hypersurfaces, 

A{V, V) + 2DvH 

is minimized hy V = iV^, • • • , F® = • • • , F") and thus to prove the main 

theorem, it suffices to verify that for alH > 0 

TT 

dtH - - nXH + — > 0. 

^ 2t- 

Fix 0 < e < T. We will apply the maximum principle to 

^ e-nXH AiL+lAPiL-fcuViiLVjiJ 

^ H H 

on the time interval [e,T). Using Lemmas 1, 5, and the inverse-concavity of the 
mean curvature, we calculate 

dtQ = dt% >AQ + 2(Vg, + 2 Q 2 + 2XH^ 

H H 

mn pQ 

^ H 

>AQ + 2(VQ,^)+2Q2. 

H 

An ODE comparison with q{t) = — 2 (t-e) which satisfies -^q{t) = 2q^[t) and 
lim q{t) = —oo shows that Q{-,t) > <z(0 ^^7 t > £■ Allowing £ —>■ 0 com- 

pletes the proof of the Harnack estimate. 


3.2. Backwards convergence. 

Lemma 6. Any weakly convex ancient solution of the mean curvature flow satisfies 

|A| < Co exp(nt) 

for t < 0 for some cq depending only on Mq. 

Proof. By the strong parabolic maximum principle, any weakly convex ancient 
solution of the mean curvature flow must become strictly convex, unless it is a 
non-moving equator. By Theorem 11 for strictly convex, ancient solutions we have 

dtH - WV.HV.H - nXH ^ 

- 

Since A = 1, we get 

dt logiL > n. 

Integrating both sides of this inequality against dt on the time interval [t, 0] gives 

H{-, t) < H{-, 0) exp(nt). 

□ 
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Lemma 7. Any weakly convex ancient solution of the mean curvature flow satisfies 

|V"i ^|2 < Cmexp( 2 ni) 

for t < 0, where Cm depends only on Mq and m. 

Proof. The proof follows by induction on m. Using the fourth evolution equation 
in Lemma 1, = A* VA and that the commutator A]T is given by 

k 

[V*^, A]T = 

j=o 

for any tensor T, we can compute the following evolution equations: 

( 1 ) 

dt\A\^ = - 2|Vy4p + 2\A\^ + 2X{2H^ - n\A\^) 

( 2 ) 

=A\V'^Af - 2|V™+^v4|2 
+ VM * V^A * V’^A * V^A 

+ AV™A * 

On the other hand, by Lemma 6 we get 

< A|Ap — 2 |VAp + Cl exp( 2 nt), 

and 

dt\yA\^ < A|VA|2 +c 2 |VA |2 +C 3 exp( 2 nt). 

Here, ci, C 2 , C 3 > 0 are independent of t. Therefore, 
dt{{t - s)|VH|2 + b\A\^) <A{{t - s)|VH|2 + b\A\^) 

+ {1 — 2b+ {t — s)c 2 )|V 2 l|^ + (ci 6 + C 3 ) exp(2nt). 

We want to apply the maximum principle on the time interval [s, s + 1]. We choose 
b large enough, independent of t, so that the coefficient of |VHp becomes negative. 
Thus the maximum principle implies that 

(t — s)|VHp < C 4 exp( 2 nt)(l — exp( 2 n(s — t))) + bc^ exp( 2 ns). 

In particular for t = s + 1 we have 

|VH|^(-,t) < C5exp(2nt). 

This verifies the bound on |VH|. Higher derivative estimates follow by induction 
and using the test function {t — s)|V"*Hp + : 

dt\V^-^Af < A\V^-^Af - 2\V^Af + Cl exp( 2 nt). 
dtlV'^Al^ < A|V’"H |2 +c 2 |V™H |2 +C3exp(2nt). 

□ 

Having established the higher derivative curvature bounds, convexity of Mt im¬ 
plies that the backwards limit of Mt is independent of subsequences. Therefore we 
have proved: 

Theorem 8. Let Mt be an ancient, embedded, weakly convex solution of the mean 
curvature flow. Then Mt converges exponentially fast in C°° to a unique equator 
M -00 as t ^ — 00 . 
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4. Classification of Ancient Solutions 

In this section, we use Theorem 8 to classify convex, embedded ancient solution 
of the mean curvature flow on 8"“''^. The proof uses the Aleksandrov reflection as 
in [4], 

We begin with some preliminaries of the Aleksandrov reflection on First, 

we will work relative to the limiting equator obtained in Theorem 8. Let E be 
an equator that bounds the open hemispheres with centers ± po, and let 
e = Opo be the unit vector in that points from the origin O to po- Let 

p{x) = (ign+i(po,x) denote the spherical distance from po to a; G The radial 

projection onto E is the map x € i-A tt{x) G i?, where tt is the nearest point 
on E to X. If X ^ ± Po, then 7r(x) is a single point. If x = ± po, then 7r(x) = E. In 
any event, given y G 7r(x), there is a unique length minimizing geodesic joining x 
to y and this geodesic must pass through ±po- 

It is convenient to make use of the ambient and define the height function 
h(x) = (x, e). The radial distance is related to the height function via 

h(x) = cos(p(x)) 

which is monotonically decreasing in p. 

Now for the Aleksandrov reflection, let V G be any unit vector that (V, e) < 
0. Let Py = V be the hyperplane through the origin with the normal vector V. 
Let = {±(x,V) > 0} denote the halfspaces with the boundary Py. For any 
subset S C write = ^TlP^. Lastly, let <5 > 0 denote the angle V makes 

with E] therefore, sin 5 = (V, —e). 

Definition 9. The Aleksandrov reflection across Py is the map defined by 

Py : X G R”+^ X - 2(x, V) V. 

This map is an (orientation reversing) isometry of R"’*’^ fixing Py and in par¬ 
ticular fixing the origin. Therefore, it induces an isometry of For x G P, we 

have (x, e) = 0 and 

h(Py(x)) = (e,x — 2(x, V) V) = 2sin5(x, V). 

In the case x G P"*", h(Py(x)) > 0, and in the case x G PnPy, h(Py(x)) = 0. 

In geodesic polar coordinates, (p,cr) G (0,7r) x P cs \{± po}, a smooth, 
closed hypersurface M C that bounds a region in Py is a smooth graph 
{f{a),a) over P if and only if its (outer) normal u satisfies (u,e) < 0 (e.q., M has 
no vertical tangents). In particular, for all e > 0 there is a ^ > 0 such that if M is 
a graph with (u, e) < —^ and N is e-close to M in C°°, then A is a graph over P. 

Theorem 10. Let Mt he a eonvex, embedded ancient solution of the mean curvature 
flow on Then Mt is a family of shrinking geodesic spheres emanating from 

an eguator at t = —oo. 

Proof. Let P = M-^o — be the limiting equator lim Mt. 

Since Mt smoothly converges to P as t t — oo, we may write Mt as the graph 
of a smooth positive function over P in the geodesic polar coordinates, Mt = 
{(/t(cr), cr) G (0, tt) X P}. We have ft —>• tt/2 smoothly and uniformly in , for any 
k, as t ^ —oo. Moreover, for 6 G (0,7r/4), Ry{E) = {{g-ao{a),a)} is a graph over 
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E. Note that Ry{E) is not a graph when S = 7r/4 (it is an equator perpendicular 
to E), so let us fix a (5o € (0,7r/4) to give us a little room. Then it follows from [4, 
Lemma 5.2] that both Mt^ and Ry{{Mt)^) are non-empty and are graphs for all 
times t G (—oo, S) with S > —oo independent of i5 S (0, i5o). Let us write gt for the 
graph of Ry{Mt) over E. 

As noted above, h(i?y(a;)) > 0 on £1+ and h(a;) = cos(/o(a;)). Thus, continuity 
implies that for e > 0 there exits an ry > 0 such that p{Ry{x)) < 7r/2 — e provided 
X G Ejj = {x € E : d{x,E OPy) > r]}. Since Mt E, we can choose < 0, 

such that d{Mt,E) < e/2 for all t < Tg; that is, p{x) > 7r/2—e/2 for all x G Mt. Now 
for X G M/" riTT”^ Erj, since Ry is an isometry, we have d{Ry{x), Ry{Tr{x))) < e/2; 
therefore, p{Ry{x)) < — e/2. Consequently, away from the strip {x G E ■. 

d{x, E n Py) < 1 ]}, we have p{Ry{{Mt)^)) < 7r/2 — e/2 and p{Mj~) > 1x12 — e/2. 
That is, Ry{{Mt)'^) > {Mt)^ away from the strip. 

Now consider the great circle C in E joining 7r(x) to Ry{7r{x)) for x G . 
From the formula 7r(x) = where x' = x — {x, e)e is orthogonal projection onto 
the hyperplane e-*- and the fact that 

we see that C lies in the plane spanned by V = V —(V, e)e and it{x). Letting y be 
the closest point to 7r(a;) in the intersection C fl we find that in fact C lies in 
the plane spanned by V and y. In particular, since 

(V,v') = cos2(5), 

the angle C makes with P is independent of the choice of x. This allows us to em¬ 
ploy the Backwards Approximate Symmetry Lemma [4, Lemma 5.1] (which applies 
whenever ft converges smoothly to C) to conclude that (possibly by decreasing Tg) 
over C 

^vii^t)y) P {Mt)y 

on the strip {x G E : d{x,Er\Py) < p}. Here, Tg depends only on \\ft — 7r/2||c2 
and therefore, Tg > —oo is independent of C. Thus we find that 

(4.1) i?v((Mt)+) > (Mt)^ 

everywhere for all t G {—oo,Tg). 

Let us now define Tg so that (—oo, Tg) is the largest interval on which the relation 

(4.1) holds. Also define T = inf Tg. We want to show that T > —oo, and hence 

se{o,So) 

that the relation (4.1) holds on the non-empty, open interval (—oo,T). To show 
T > —oo, we apply the maximum principle: Recall that both Mt^ and Ry{{Mt)^) 
are non-empty and are graphs for all t G (—oo, S) with S > —oo independent of 
S. Since the relation (4.1) is true on (—oo, min{T 5 , S'}), the maximum principle 
applied in the time interval [minlT^, S}/2, S) ensures that 

Rymt)p > (Mt)^ 

for all t G (—oo, S) and any J g (0, Jq)- Hence T > S. 

To complete the proof, we use [4, Proposition 5.3] (which applies in any dimen¬ 
sion) to conclude that Mt is a geodesic sphere for all t G (—oo,T) and thus for all 
negative times by the uniqueness of solutions. □ 
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